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The idea of Dark Matter (DM) with self interaction was invoked to resolve a 
number of discrepancies between the simulation based predictions by collisionless 
cold DM and the astrophysical observations on galactic and subgalactic scales. Ev¬ 
idences for self interaction would have striking implications for particle nature of 
DM. In order to reconcile such astrophysical observations for self interaction with 
particle properties for DM, we consider the general scenario of self interacting Dirac 
fermionic DM, y. Also since the exact particle physics model for DM is yet to be 
probed, we simply adopt the effective model independent framework for DM self 
interaction which occurs via the most general effective 4-fermion operators invariant 
under both Lorentz and CPT transformations. Prom the thorough investigation of 
the interrelations among the parameters in this framework, namely, the effective DM 
self couplings (Gj), DM mass (m^) and relative velocity (urei), it can be inferred that 
Gi decrease with increasing for a given DM self interaction strength. Moreover, 
for few types of effective operators the values of Gi fall off with increasing Urei while 
they remain roughly constant for a wide range of Urei for other cases. In addition, 
the parameter space in this framework is constrained by the claimed observational 
results of ^ on cluster scales (Abell 3827, Bullet Cluster) after averaging the DM 
self interaction cross sections over DM velocity distribution in the cluster. This 
puts interesting constraints on the values of effective DM self couplings for different 
fermionic DM masses for various effective operators (scalar, vector, etc.) of DM self 
interactions in this scenario. Some other implications of DM effective self interaction 
are also discussed in this model independent framework. 
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I. INTRODUCTION 


It is now very well established from various astrophysical and cosmological evidences M 
that a large fraction of the universe’s mass consists of unknown non-luminous matter, namely, 
Dark Matter (DM). The Planck satellite-borne experiment jd] estimates from the anisotropy 
study of the expected isotropic nature of Cosmic Microwave Background Radiation (CMBR) 
that ~ 26.8% of the total energy content of the universe is dominated by DM. Since there is 
no candidate of DM in the Standard Model (SM) of particle physics, it is believed that there 
exists a new physics sector beyond the SM. Although the nature of DM remains unknown still 
now, it is conjectured that the weakly interacting massive particles (WIMPs) being stable, 
non-relativistic, electrically neutral and colourless are very well-motivated candidates of DM. 
The mass of WIMPs can vary from a few GeV to TeV with suitable strength of interaction to 
yield observed DM relic abundance depending on various models of particle physics 

On the other hand, the formation and evolution of structures in the universe can also 
predict precise nature of DM sector. The large scale structure of the universe can almost 
be successfully accounted for in case of cohisionless cold dark matter (CCDM) paradigm. 
But such scenario fails in explaining the small-scale structure of the universe. There are 
persistent problems such as core-vs-cusp problem 20|, missing satellite problem 2l(], too- 
big-to-fail problem 2^ etc. which arise in the CCDM scenario while predicting the precise 
nature of the small scale structures in the universe. The central density prohles of the DM 
halos surrounding the dwarf galaxies indicate flat cores from the observation of galactic 


19|. 


rotation curves 


steep cusps [28 


23 

30 


of galaxy clusters 


27l | while numerical simulations considering CCDM model predict the 


31 


Similar discrepancy has also been noted for the case of observations 


33|. This discrepancy known as the core-vs-cusp problem cannot be 


resolved by effects such as supernova feedback etc. Another problem, namely, the missing 


satellite problem is the discrepancy of ~ 0(10) magnitude be 


observed numbers of satellites within the Milky Way 


£, 34 


ween the expected and the 


Although recently the 


Sloan Digital Sky Survey has discovered many faint galaxies with much precision, there 


has been a multiplying factor o 
and undiscovered dwarf galaxies 
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5 — 20 discrepancy between the numbers of known 


39|]. Therefore there exists an unresolved issue on the 


predicted number of Milky Way subhalos in CCDM scenario. Moreover, there exists another 
interesting discrepancy whereby the estimations of numerical values of host subhalo mass 














































3 


and density become so big that they fail to confront the observational evidences, known 
as the too-big-to-fail problem. The largest velocity dispersions of the brightest Milky Way 
dwarf spheroidal (dSph) galaxies observed from the rotation curves of these galaxies may 
indicate that these galaxies are hosted by the largest subhalos in the Milky Way halo. But 
the central densities of such highly massive host subhalos as predicted by the CCDM-only 
simulations are too large and also the masses of such subhalos become too huge to host 


those above-mentioned brightest dwarf satellite galaxies 


22 


40 


4l| . The maximum circular 


velocities of these satellite galaxies (or dSphs) are observed to be much less lower than those 
estimated by numerical simulations in CCDM paradigm. This too-big-to-fail problem may 
be attributed to have similar resolution like that of the core-vs-cusp problem since such 
massive subhalos can be made consistent with the observed dSphs if the central densities 
are reduced. 

With these long-established discrepancies involved in the CCDM paradigm, a suitable 
modihcation of such scenario is much required. Inclusion of warm DM j^, i.e., DM with 
some kinetic energy was the hrst attempt towards such modihcation to resolve the missing 


satellite problem. However, the other severe prob 
warm DM 3, 44|. There are several valid attempts 


ems cannot be fully resolved with this 


45 


52| to alleviate such problems. T 


include the proposition of self-interacting DM as a solution to all such discrepancies 


. T1 

lese 

45 I, 

46|. 


The energy transfer from the outer hotter region in the halo to the central colder region can 
form a core structure within a central region of the halo in self-interacting DM scenario. 
Also, in this case the number of Mil^ Way satellite galaxies is much reduced due to the 
collisional stripping of the subhalos [53]. 

Since the beginning of the proposition of the self-interacting DM, various numerical simu¬ 


lations 


^58 1 have tried to estimate its value for different astrophysical observations. Early 
simulation had predicted ^ ~ 1 — 10 cm^/g required to yield the observed flattening of the 
central densities of dwarf galaxies SJ]- Moreover, the abundance of Milky Way subhalos is 
considerably reduced for ^ ~ 10 cm^/g. Other simulation 59|] had found that the value of 
^ to be ^ 0.1 cm^/g to avoid core formation in halos of galaxy clusters, which violates the 
gravitational lensing observational constraints for the cluster CL 0024-1-1654. The param¬ 
eter space of ^ ~ 0.3 — lO'^ cm^/g is ruled out requiring the non-evaporation constraints 
of elliptical galaxy halos within hot cluster halos [GOj whereas the cluster ellipticity lim¬ 
its require — ^ 0.02 cm^/g {^1- More recent numerical study on clusters, however, have 
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found no significant evaporation of subhalos for — ~ 1 cm^/g 5a, 1^- The current simula¬ 


tions 


55 


57l | hint the upper and lower limits of self-interaction as — ~ 0.1 — 10 cm^/g to 


resolve the core-vs-cusp and too-big-to-fail problems on dwarf scales. Although there are 


controversies 


regarding such predicted lower limit of self-interaction 


0.1 cm^/g). 


On the other hand, the observational constraints from cluster and Milky Way scales restrict 
the self-interaction as — ^ 0.1 — 1 cm^/g 

m ^ ® 


55 


m- 


In a very recent analysis 62|], it is claimed from the observation of a particular galaxy in 
the cluster Abell 3827 with the Hubble space telescope imaging that an interesting offset of 
~ kpc (68% CL) between its DM and its stars (normal baryonic matter) is observed. 

This offset is thought to be a consequence of self-interaction of DM (^ ~ (1.7 ± 0.7) x lO”"^ 
cm^/g), and a lower bound on ^ ^ 10“^ cm^/g is given from their analysis 62|. On 
the other hand, using a very different kinematic analysis for the same galaxy, the authors of 
Ref. 63j have obtained another lower bound on ^ ^ 1.5 cm^/g (more precisely its order 
lies ~ 1.5—3.0 cm^/g). The analysis of the other astrophysical bodies such as “bullet cluster” 
(IE 0657-56) [6^ sets the limit (upper) on the quantity ^ to be less than ~ 1.25 cm^/g (68% 
CL) considering the effect of drag force arising from DM particle collisions on the subcluster 
halo (particle momentum exchange is isotropic). A more stringent limit, ^ < 0.7 cm^/g 
(68% CL), is also set by considering the effect of mass-to-light ratio of the subcluster due to 
the scattering of DM particles (particle momentum exchange is directional) 6^. Another 
tight constraint on self-interaction has been derived from the weak lensing analysis of the 


observational data of the giant bullet cluster collision 1E0657-558 65| . Also the studies 
on 72 clusters collisions data put an upper limit on —, — < 0.47 cm^/g with 95% CL [66|. 


Clearly the obtained result of Ref. 


66 | is in tension with the previous obtained bound from 


the observation of “bullet cluster” (IE 0657-56). 

Although various aspects of self-interacting DM scenario have widely been investigated 
in numerical simulation as well as in astrophysics, those for particle physics ^ are com¬ 
paratively less available in literature. Also, studies on both DM velocity dependent and 
independent self-interactions have also been pursued. While early studies mainly focus 


on the velocity-independent D 


velocity-dependent scenario 68 


VI self interaction, there have been recent developments in 
which can evade observed bounds by yielding different 


DM self scattering cross sections on various scales (dwarf galaxy, cluster etc.) and are thus 
interesting to study. There are a variety of other proposed models for DM self interactions: 
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mirror DM 


76 


78| . atomic DM [79|, [8^, DM coupled to dark photon 0 which is originated 


hidden symmetry 


68 


7 ^. 


Although there are a plethora of particle physics models for WIMP which can explain 
various astrophysical as well as experimental observations, however, there are still lack of 
experimental or observational evidences to precisely distinguish the correct particle physics 
model for DM physics. It may be possible to detect DM in future with direct or indirect 
detection techniques but still these early observations may only provide the information 
about the properties of DM particles but not the underlying theory. Therefore in order 


to avoid such theoretical bias, model indepenc 


particularly interestin g in rec ent studies 
There are studies 


81 


86 


ent studies of DM phenomenology become 


87 


97| from the aspects of such model independent ways on DM 


phenomenology to relate the observed DM relic density, direct & indirect detections and 
collider signatures. In this paper we also consider a model independent WIMPs rather than 
choosing any particular model. We follow the approach of effective 4-fermion interactions 
among the WIMPs for various combinations of spins and interaction forms 84, l96|]. For this 
study we have simply considered a Dirac fermionic WIMP and construct the general effective 
4-fermion operators for DM-DM self-interactions. A detailed investigation and analysis of 
parameter space is performed thereafter. Moreover, we also evaluate and interpret the 
constraints for the self-interaction of WIMPs for each case of spin and interaction forms 
from cluster observational results. 

The paper is organised as follows. In Sec. Hi] the effective Lagrangians for DM self¬ 
interaction in this framework are briefly discussed. The calculational results for effective 
DM self-interaction cross sections are mentioned in Sec. urn The following section (Sec. HV]) 
is devoted to a detailed discussion on the parameters of this effective DM self-interaction 
scenario. The effective self-couplings are constrained for a wide range of DM masses from 
several observational bounds on galaxy cluster scales in Sec. jVj In Sec. jVT] we hnally sum¬ 
marise our work and some interesting conclusions are drawn. 


II. EFFECTIVE LAGRANGIANS FOR DARK MATTER SELF-INTERACTION 

We consider the case of a single Dirac fermionic WIMP, % for pursuing this analysis. In 
order to begin our study we take effective operators between DM particles without consid- 
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ering any particular preferred model for DM. However, we make some assumptions at this 
point. First of all, the WIMP is assumed to be the only new particle beyond SM at the 
electroweak scale. This implies that the calculations performed hereafter remain unaffected 
by possible other effects (resonances etc.) and self interactions can thus safely be described 
by effective held theoretical framework. Secondly we also assume, for simplicity, that the 
considered 4-fermion interactions among the DM fermions is governed by only one form of 
interaction (scalar, vector, tensor, etc.) among the possible set of effective operators. There¬ 
fore the possible effective Lagrangians involving effective self interaction operators among 
four such WIMPs (x) can be given by. 


Scalar (S) ; 

= Gs XXXX , 

(ii.i) 

Pseudoscalar (P) : 

>Cp = Gpxl^xXlbX , 

(11.2) 

Vector (V) : 

£v = Gv XYXX1^^X , 

(11.3) 

Axialvector (A) : 

= Ga Xl^l^XXlpGhX , 

(11.4) 

Tensor (T) : 

£t = GtXC^^''XW^ ivX , 

(11.5) 

Scalar-pseudoscalar (SP) : 

>Csp = Gsp XXXilbX , 

(11.6) 

Pseudoscalar-scalar (PS) ; 

£ps = Gps X^lbXXX ) 

(11.7) 

Vector-axialvector (VA) ; 

£va = Gva X7^XX7m75A , 

(11.8) 

Axialvector-vector (AV) : 

Cm = Gav X7^75XA7mX > 

(11.9) 

Alternative tensor (T) : 

Cf = Gf e^''P''x<yt^uXX(^paX , 

(ILIO) 

Left handed-left handed (LL) ; 

^LL = Gll X7^(1 - 75)xX7m( 1 “ 75)x , 

(11.11) 

Right handed-right handed (RR) ; 

^RR = Grr x7^(1 4- 75)xX7m( 1 + 75)x , 

(11.12) 

Left handed-right handed (LR) : 

Clk = GlrX1^{1 - 75)xX7m( 1 + 75)x , 

(11.13) 

Right handed-left handed (RL) ; 

^RL = Grl Xl^il + 75)xX7m(1 - 75)7 , 

(11.14) 


where the effective coupling constants Gs are the real number with mass dimension —2. In 
the above 7 ^, 75 are the usual gamma matrices and are the antisymmetric rank -2 

tensor and Levi-Civita symbol respectively. We tag the above interaction types based on 
the 2 —)■ 2 self-scattering processes of fermionic DM particles, namely, XX XX using 
effective theoretical approach. It should be noted that the considered effective operators 
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Scalar 

Pseudoscalar 

V ector 
Axialvector 
T ensor 
Chiral 

Alternative Tensor 
Scalar — Pseudoscalar 

V ector — Axialvector 


FIG. 1. Feynman diagram for self interaction of fermionic dark matter (x) in effective theory. 
See text for more details. 


shown in Eqs. Ill.ltlll.l^ are invariant nnder both CPT and Lorentz transformations. Also, 
the hrst hve operators in Eqs. III.ltlII.5l are separately invariant under C, P and T trans¬ 
formations. In addition, they also respect the Hermiticity requirements. All of the above 
listed operators do not guarantee independent self interaction cross sections. For example, 
the effective Lagrangian for scalar-pseudoscalar (SP) interaction, £sp of Eq. III.61 and that 
for pseudoscalar-scalar (PS) interaction, £ps of Eq. III.71 are similar in nature for the case 
of effective 4-point diagrams of DM self-interaction and can yield similar results. Simi¬ 
larly, the effective Lagrangian for vector-axialvector (VA) operator of Eq. III.81 and that for 
axialvector-vector operator of Eq. III.91 result equal squared transition amplitude (or squared 
scattering amplitude) for effective self-interacting diagrams. Moreover, the last four effective 
Lagrangians of Eqs. Hl.llllII.l^lII.ldl and lll.ldl namely, £ll, -Crr, £lr and £rl respectively, 
represent chiral (C) interactions among DM particles and also produce squared amplitude 
terms for DM self-scattering diagrams similar to each other. Hence from the above listed 
set of operators we are left with nine operators that can result in producing independent 
self scattering cross sections in case of Dirac fermionic DM. Note that the chiral interaction 
operators can be derived from suitable combinations of other listed operators. Also there 
may be some other forms for alternative tensor (T) in Eq. III.101 such as 
which are actually equivalent to xa^yXX^^puX- 

The basic Feynman diagram for the aforementioned effective operators is depicted in 

Fig. [H 






















III. EFFECTIVE DARK MATTER SELF INTERACTION CROSS SECTIONS 


In order to continue our analysis we calculate the DM self scattering cross-sections for 
self interaction considering the above mentioned operators and get the following 
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(IIM) 

(111.2) 

(111.3) 

(111.4) 

(111.5) 

(111.6) 

(111.7) 

(111.8) 

(111.9) 


where s is the usual Mandelstem variable, my_ is the mass of fermionic WIMP y. Now 
since in the non-relativistic limit the relative velocity can be written in terms of mass of the 
WIMP (x) and the Mandelstem variable s as, 

^s{s - Am^) 


(s — 2m^) 

Hence in the lab frame s can be rewritten in terms of the relative velocity, Prei as 

1 


(III.IO) 


s = 2m^ 1 + 


y/T^V. 


2 

rel, 


(111,11) 


Also in this context it should be mentioned that if the velocity of two colliding DM particles 
are Vi and V 2 , then Prei can be written as 


Vrel = \Vl - V2\, 


( 111 . 12 ) 
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in the non-relativistic limit, while in the case of relativistic velocities, the relative velocity 


between such DM particles is given by 


'^^rel 


98| 


(d - V 2 ) 


2 {viXV 2 )'^ 

^2 


(III.13) 


Velocity of light c is chosen to be unity and hence Vrei is dimensionless in our later calcula¬ 
tions. 

Substituting the expression of s of Eq. IIII.lll Eqs. Illl.ltllll.^ can be recasted as 
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CTt, self — 2?Tr^ 


3 + 
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-^rel 


1 + 
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Cri^,self = CT^P self — (^PS,se\i — 2m. 

lOTT 
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(111.14) 

(111.15) 

(111.16) 

(111.17) 

(111.18) 

(111.19) 
, (III.20) 

(111.21) 

(111.22) 


In the above, since asp^seii and aps^seii yield similar results, these two terms are denoted 
by cr;^,seif (say). Similarly, UM.seif represents the terms uy^^seif and crAy,seif while ac,seif is 
for the chiral (LL, RR, LR, RL) ones. Now, it can be noted from Eqs. IIII. 1411111.221 that 
the DM self scattering cross sections depend only on three parameters, namely, the effective 
couplings, DM relative velocity and DM mass. 
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IV. PARAMETRIC DEPENDENCE OF DARK MATTER EFFECTIVE SELF 

INTERACTION CROSS SECTIONS 

In order to inspect the general nature of Eqs. IIII.14llnL22l all of the above equations 
(Eqs. IIII.141liri.22l) can be recasted in terms of various types of scattering cross-sections (cXi) 
in a general form as, 

a, = QG^,2mlf,{v,,,), (IV.l) 

where GjS are the effective couplings and /i(nrei) are different functions of DM relative veloc¬ 
ity for different operator types. In the above Ci are constants which depend on phase-space 
factors for different types of interactions mentioned before. Now it can be seen from Eq. IIV.II 
that for any type of effective self-interacting DM scenario, the self-scattering cross-section is 
proportional to the square of both effective self-coupling and DM mass, i.e., Gfrn^. Thus for 
any type of self-interaction in such effective interaction theory, a oc G^m^ for a particular 
DM relative velocity. Since from various astrophysical as well as cosmological observations 
one can put limits on the factor, this could restrict the limit on the parameter space 

containing effective self-coupling, DM mass and its relative velocity. In this context it is 
worth important to note that the DM relative velocity depends on the thermal distribution 
of DM particles. If the DM relative velocity is estimated from such thermal distribution of 
DM, one can, in principle, estimate the limits on the effective self-coupling for different DM 
masses from the precisely measured DM self-scattering cross-section. 

In order to elaborate how the self-interaction parameters depend on each other in this 
effective theory, we make a 3-D colour coded plot in Fig. [2] where the variations of ef¬ 
fective self couplings {Gi) with the simultaneous variations of other parameters, namely, 
dark matter relative velocity (vrei) and dark matter mass (m^) are furnished. In Fig. [2] 
we have estimated such variations for the above-mentioned nine independent operators of 
self-interaction processes, namely, for scalar (S), vector (V), pseudoscalar (P), axialvector 
(A), scalar-pseudoscalar (SP), vector-axialvector (VA), tensor (T), alternative tensor (T) 
and chiral (C) interactions. The colour index in Fig. [2] denotes the DM mass where the DM 
mass varies from bluish coloured region to yellowish region as it increases. For all of the 
plots in Fig. [2] the mass of DM particles has been varied from 10 MeV to 1 TeV as shown 
by the colour coded index. However, we can, in principle, choose any other mass ranges 
of DM but the resulting plots will be similar in nature. Moreover, we have considered two 
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FIG. 2. Variations of effective self couplings (G,) with both dark matter relative velocity (vrei) 
and dark matter mass (m^) for different possible types of self-interactions discussed in the text. 
The colour gradient indicates the variation of dark matter mass. The blueish coloured region in 
the colour bar denotes the lower mass zone while the yellowish one is for the higher mass. The 
two band-like regions in each plot are for two considered values of self-interaction cross-sections, 
namely, for ^ = 5 x 10^ GeV~^ (10.0 cmf/g) and 5 x 10“^ GeV~^ (l.Q x 10“® cmf/g). See text 
for more details. 


different values of self-interaction cross-sections, namely, ^ = 5 x 10^ GeV ^ (10 cm^/g) and 
5 x 10“^ GeV“^ (1.0 x 10“^ cm^/g) for each of these plots. 0 This is because of the fact that 


^ Note that 1 GeV ^ = 2 x 10 ^ cm^/g, 1 barn/GeV = 0.6 cm^/g 
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the self-interaction is not precisely measured since and the measurements of the factor, ^ 
from various observations vary much significantly (few orders of magnitude) due to various 
astrophysical uncertainties. Hence we restrict the values of the factor, ^ to those which 
may be the predicted extreme limits on self interaction from various studies 67(]. We can, in 
principle, choose any other limiting values of ^ which were estimated to confront several 
astrophysical observations such as bullet cluster etc. This would rather change the limits on 
the effective self-couplings. But the main motivation for this plot remains unchanged. The 
two band-like structures in each frame of Fig. [2] are due to such two choices of In each 
plot the upper band corresponds to the chosen higher value of ^ (10.0 cm^/g or 5 x 10^ 
GeV“^) while the band on the lower side is for the lower value (1.0 x 10“® cm^/g or 5 x 10“^ 
GeV“^). In Fig. [2] we see that for a particular DM mass the effective self-couplings remain 
almost constant with a broad range of DM relative velocities for the self-couplings to be 
vector (V), axialvector (A), pseudoscalar (P), tensor (T) and chiral (G) in nature. On the 
other hand, for the cases of effective self-couplings due to scalar (S), scalar-pseudoscalar 
(SP), vector-axialvector (VA) and alternative tensor (T) interactions, these effective self¬ 
couplings decrease with increasing DM relative velocities as shown in Fig. [2l However, at 
a very large DM velocities, i.e., for the case of relativistic self-collisions of DM particles, 
the effective self-couplings decrease very steeply with DM relative velocity for all types of 
self-couplings. It is clear from Fig. [2] that these above-mentioned variations of different 
self-couplings with DM relative velocity are similar for any chosen values of DM mass and 
self-interaction cross-section. It is also obvious from this hgure that for a particular values of 
?;rei and if the DM mass becomes very high (towards yellowish coloured zone in Fig. |2]), 
the value of self-coupling decreases compared to the case where DM mass is smaller (towards 
bluish coloured zone in Fig. [2]). 


We hnally estimate the dependences of various types of effective self-couplings (Gj) on 
different mass scales of fermionic DM (m^). In order to compute such dependencies we 
again only consider the two recently reported values of DM self-interaction on Abell 3827 
observation. Needless to mention here that we can, in principle, consider any reported 
DM self-interaction values from astrophysical or cosmological observations as well as from 
numerical simulations other than those chosen in this work. But this may change the 
constraints on the effective self-couplings from that obtained in this work. Since there 
are conflicts among different reported self-interaction values, we keep our options open in 
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Dark Matter Mass (GeV) 

FIG. 3. Plot showing the allowed parameter spaee in the plane containing effective self couplings 
(Gi) and fermionic dark matter mass (m^) for different possible self interaction types. The black 
dotted line in each plot denotes dark matter relative veloeity, Uj-ei = 0.999c (relativistie collision) 
while the solid black line is for Vrei = 0 (non-relativistie collision). The red and magenta eoloured 
regions in each plot is for two different self interaetion eross section values, namely, for ^ = 5 x 10^ 
GeV~^ (^10.0 crnf/g) and 5 x 10“° GeV~^ (l.t) x 10“^ cmf/g) respectively. The green and blue 
region denote the allowed and excluded parameter space for the above two values of ^ whereas the 
yellow region is allowed for ^ = 0.05 GeV~^ but excluded for the other ^ value. See text for more 
details. 

this work. We show the allowed zone of all possible types of effective self couplings (Gj) 
for different dark matter masses {m^) in Fig. 121 We would like to mention at this point 
that similar to Fig. [2l here also the considered interactions include scalar (S), vector (V), 
pseudoscalar (P), axialvector (A), scalar-pseudoscalar (SP), vector-axialvector (VA), tensor 
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(T), alternative tensor (T) and chiral (C) types. In all of these plots we have taken into 
consideration the two extreme values of DM relative velocity (urei), namely, Vrei = 0.999c 
and 0 (c is velocity of light in vacuum). These velocities correspond to the relativistic and 
non-relativistic regimes of DM self-scattering respectively. In Fig. [3] the black dotted lines 
indicate dark matter relative velocity, Vrei = 0.999c while the black solid lines are for the 
other one, Vrei = 0. The red and magenta coloured regions in each plot of Fig. |3] denote the 
allowed zone for two considered self-interaction cross section values, namely, for ^ = 5 x 10^ 
GeV“^ and 5 x 10“^ GeV“^ respectively. Also in each plot of Fig. [3] the green and blue 
regions denote the parameter space for values of ^ < 1.0 x 10“^ cm^/g (5 x 10“^ GeV“^) 
and ^ > 10 cm^/g (5 x 10^ GeV“^) respectively, whereas the yellow region is allowed for 
the values of DM self-interaction in between the chosen ones, 1.0 x 10“^ cm^/g (5 x 10“^ 
GeV“^) < ^ < 10 cm^/g (5 x 10"^ GeV“^). We see in Fig. |3] that the allowed zone of 
effective self couplings, Gi decrease with increasing DM mass. This simply corresponds to 
the fact that self-interaction of more massive fermionic DM particles happens through a 
heavier messengerparticle while a lighter messenger is responsible for the self-interaction 
of less massive fermionic DM. The black solid and dotted lines in each plot of Fig [3] gives 
the upper and lower bounds on Gi for different DM masses, for a particular value of 
self-interaction. Since we have considered two values of — for our calculations, there should 
be two black solid and dotted lines in each plots of Fig. [3l But in Fig. |3] it can be seen 
that there in no black solid line (urei = 0) appearing in few frames, namely, for the cases 
of scalar (S), scalar-pseudoscalar (SP), vector-axialvector (VA) and alternative tensor (T) 
interactions. So, there is no restriction on upper bounds for those interactions. This is 
due to the nature of calculated theoretical forms of /j(nrei) in self-annihilation cross sections 
for those types of interactions. Therefore, the constraints on such types of self-interaction 
is relaxed to higher values if DM particles interact much more non-relativistically. Hence 
the blue regions are shifted to much higher values of Gi and do not appear in those plots. 
Also since the magenta and red regions get partially overlapped with each other due to such 
shifts of black solid lines in those plots, the yellow regions extend up to large values and get 
overlapped by these magenta and red regions. Thus only lower bounds on Gi for different 
exist in those cases, which are denoted by black dotted lines. 
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V. CONSTRAINING EFFECTIVE DARK MATTER SELF COUPLINGS FROM 

OBSERVATIONS IN CLUSTER SCALES 


In the realistic cases, the DM particles do not traverse in space with a hxed velocity but 
follow a velocity distribution such as Maxwell-Boltzmann velocity distribution. Therefore 
the DM scattering probability in the DM halo is not hxed by only a single value of DM rela¬ 
tive velocity, t’rei but by the convolution over different DM densities and velocities, precisely 
determined by N-body simulations. Here for simplicity we make an average of the scatter¬ 
ing cross section over the velocity distributions of interacting DM particles. This velocity- 
averaged self-interacting cross section of DM ((a)) can be approximately treated as the 
measured quantity on which various astrophysical observations put different constraints 75 1. 
Since the WIMPs are assumed to be thermally distributed we consider Maxwellian velocity 
distribution for the DM particles of the form 


Hv)^ 



(V.l) 


where Vq is the most probable speed of DM. The choice of Vq depends on the different 
characteristic velocity scales of different halo sizes. If the initial DM velocities are vi and 
V 2 , this averaging over DM velocity distributions of such DM particles gives. 


(cr) = / (fvid^V2J^{vi)J^{v2)(T{\vi-V2\) 


(Pvid?V2 _„ 2 / 2 _^ 2 ; 2 ^ 

(7ra^)3 

(27rr'Q)3/2 


e CT(lftel) , 


(V.2) 


where Urei = l^i — frl is the DM relative velocity as discussed earlier. The characteristic 
DM velocity scale in galaxy clusters is approximately ~ 1000 km/s whereas for Milky Way 
galaxy such characteristic scale becomes ~ 100 km/s. On the other hand, in case of dwarf 
galaxy the characteristic velocity takes the value of ~ 10 km/s, comparatively much smaller 
than the aforementioned velocity scales of other systems. We choose the value of Vrei to 
be ~ 1000 km/s or ~ 3.3 x 10“^ (in the unit of velocity of light in vacuum, c ~ 3 x 10^ 
km/s) since we consider the reported observational bounds on DM self-interaction in cluster 
scales. To proceed we compute the velocity-averaging of different DM self-interaction cross 
sections given in Eqs. IIII.141lni.22l The computed values of velocity-averaging of functions 
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such as , ^ ^ and ^ ^ ^ appearing multiple times in different forms of /i(nrei), come 

out to be ~ 1.00002 and ~ 0.5 respectively. In this method the values of (a) for different 
types of aforementioned self interactions are derived for a typical cluster velocity scale of 
~ 1000 km/s. The computed ratio of (a) to DM mass (■^) in this effective framework 
may be safely treated as the observed constrained quantity for self interaction in 

cluster scales. Therefore the calculated quantity ^ may be used to put constraints on the 
parameter space of this framework containing effective self-couplings (Gj) and DM mass 
(m^). 

Needless to mention at this point that the above-mentioned velocity-averaging is not 
relevant if the cross section is independent of DM velocity whereas this averaging becomes 
more important in case of strongly velocity-dependent cross sections (for example, resonance 
effects etc.). 


In order to give a much more clear picture of the constrained parameter space, we have 
calculated the effective self couplings constrained by various observations of galaxy clusters 
for different DM mass ranges and shown the constraint parameter space containing the 
effective self couplings (G*) and the DM mass {m^) in Fig. IH Clearly, each frame of Fig. |4] 
represents a particular kind of DM self-interaction scenario discussed earlier. In each plot 
of Fig. |4]the solid green line represents — = 1.7 x 10“^ cm^/g from the recent observation 
(central value) of Abell 3827 cluster 62j while the solid red line denotes for the another 
claimed limit on namely, for ^ = 1.5 cm^/g from the study of the same cluster 63|. On 
the other hand, the upper limit on DM self interaction (^ < 0.47 cm^/g) from the analysis 
of observational data of 72 galaxy clusters 66(] is shown by the solid blue line in each 
plot whereas the black solid line is for the conservative upper limit on DM self interaction 
(^ < 1.25 cm^/g) from the observation of bullet cluster (IE 0657-56) 6^. The total allowed 
region of parameter space for ^ < 0{1) cm^/g as reported by various cluster observations 
is shown by the yellow coloured region in each plot of Fig. 01 Since different observational 
as well as simulation-based results in the recent past for the cluster scales impose severe 
constraints on the quantity ^ so that its numerical value cannot exceed ~ 0(1) cm^/g, the 
parameter space outside these yellow regions in Fig. 0] is also not viable for those results. 
This also implies that for any limit of ^ to be less than 0{1) cm^/g, the projected values of 
effective self couplings lie on the yellow regimes for the range of DM masses shown in this 
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10 "'' 10 ° 10 ^ 10 ^ 10 ° 10 "'' 10 ° 10 ^ 10 ° 10 ° 10 "'' 10 ° 10 ^ 10 ° 10 ° 


Dark Matter Mass (GeV) 

FIG. 4. Constraining various possible types of effective self couplings (Gi) for different fermionic 
dark matter masses (my') from the observational results of galaxy clusters. In each plot the solid 
green line denotes the reported limit on the the ratio beween dark matter self-interaction and mass, 
^ = 1.7 X 10“^ cmf /g from the recent observational data of Abell 3827 cluster while the solid 
red line is for another claimed limit on namely, ^ = 1.5 cmf/g from the observation of same 
cluster. On the other hand, the upper limit on dark matter self interaction (^ < 0.47 cm)/g) from 
the recent study of 72 galaxy clusters is shown by the solid blue line in each plot whereas the black 
solid line denotes the conservative limit (upper) on such self interaction (^ < 1.25 cm?/g) from 
the bullet cluster IE 0657-56 observation. The yellow region represents the total allowed region of 
parameter space for ^ < 0(1) err?/g reported by various cluster observations. See text for more 
details. 


figure. Therefore any projected limit on DM self-interaction for cluster scale observations 
can, in principle, be translated to the the constrained region of parameter space containing 
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Gi and It is obvious from Fig. 0] that for a given DM mass the variations of the 

constrained couplings Gi are within (9(10) among different types of couplings except for the 
cases of scalar-pseudoscalar (SP) and alternative tensor (T) interactions where the calculated 
bounds are found to be somewhat relaxed to much higher values of couplings. It can be 
inferred from Fig. IHthat the obtained values of DM effective self couplings constrained by 
the cluster observations for each considered DM mass scale are much higher by a few orders 
of magnitude than those of DM effective couplings with SM fermions ?? required by Planck 
relic density constraints. 


VI. SUMMARY AND CONCLUSION 

The concept of DM self interaction was hrst proposed to alleviate myriads of astrophys- 
ical problems. It also plays very signihcant role for the formation and growth of structures 
of the dwarf galaxies as well as of the galaxy clusters. In this work, we discuss the possibil¬ 
ity to constrain DM self interaction using the effective interactions among self interacting 
DM particles. For this we have chosen a general Dirac fermionic DM and considered all 
possible effective 4-fermion interaction operators among such DM particles. We have found 
nine possible types of such effective self interactions which yield independent forms for 
DM self scattering cross sections. These interactions include scalar (S), pseudoscalar (P), 
scalar-pseudoscalar (SP), vector (V), axialvector (A), vector-axialvector (VA), tensor (T), 
alternative tensor (T) and chiral (C) types. Using these effective interactions for the cal¬ 
culation of DM self interaction cross sections, it is found DM that for each type of those 
interactions, only three parameters, namely, DM effective self coupling (Gj), DM mass {rriy^ 
and relative velocity {vrei) are required. A detailed study on how these parameters depend 
on each other has been performed. From such study it has been shown that for a given value 
of DM mass and self interaction, the DM effective self couplings (Gj,^ = V,A,P,T, C) do 
not change considerably for a broad range of DM relative velocity in the cases of vector, 
axialvector, pseudoscalar, tensor and chiral types of interaction. On the other hand, for the 
cases of scalar, scalar-pseudoscalar, vector-axialvector and alternative tensor interactions, 
the DM effective self couplings (Gj,i = S,SP,VA,T) decrease with the increment of Urei- 
However, for any type of interaction there is a rapid decrement of Gj with increasing Urei 
near the relativistic regime of DM self interaction. This happens possibly be due to the 


19 


strong enhancement of self interaction cross section which, in turn, results a reduction of 
self coupling near such very high relative velocity. Likewise, the self couplings, Gi for a given 
value of ^ fall off with increasing DM mass (m^). These decrements of Gi with increas¬ 
ing follow a power law behaviour with index of power law being —0.5. We have also 
shown the allowed values of Gi for a range of for two extreme cases of Urei, namely, for 
^rei = 0 (non-relativistic collision) and 0.999c (relativistic collision). From such study it can 
be inferred that the constraints on Gi for different become tighter when DM particles 
collide (or interact) with each other relativistically while such constraints are relaxed for 
non-relativistic interaction between DM particles. However, the scale of such relaxations 
precisely depends on the nature of self couplings. The relaxation is comparatively low for 
tensor interaction among all interactions whereas for the cases of scalar, scalar-pseudoscalar, 
vector-axialvector and alternative tensor interactions it becomes so many orders of magni¬ 
tude large that there is no preferred limit for such non-relativistic interaction. In addition, 
it is to be mentioned in this context that these constrained values of Gi shift towards higher 
values if DM self interaction becomes stronger (i.e., self-scattering cross section is higher) 
keeping other parameters constant in this effective model and it also follows conversely. 

In order to get parameter space consistent with various astrophysical bounds for self 
interaction, we require the characteristic velocities for different size halos which may, in 
turn, determine the DM velocity distributions. Therefore, in order to derive the limits 
on Gi for different from the observed results of galaxy clusters, we need to perform 
velocity-averaging of the obtained velocity dependent self scattering cross sections over a 
Maxwellian velocity distribution with most probable speed being the characteristic cluster 
velocity, vo ~ 10^ km/s. Finally few observed constraints on cluster scales are put into 
such resulting velocity-averaged cross section terms ((cr)). We have found that for the 
considered choice of being 10 MeV to 1 TeV, the values of Gi,i = V, S, P, A, VA, T, C 
roughly he between ~ 10^ GeV“^ to ~ 1 GeV“^ and those for Gi,i = SP,T become ~ 10^ 
GeV“^ to ~ 10^ GeV“^ consistent with — = 1.5 cm^/g from Abell 3827 cluster. On the 
other hand, for the requirement of ^ = 1.7 x 10“^ cm^/g from another analysis, it requires 
Gi, i = V,S, P, A, VA, T, G to lie between ~ 10 GeV“^ to ~ 0.1 GeV“^ and Gi, i = SP, T to 
take values between ~ 10^ GeV“^ and ~ 10 GeV“^. Although the results are shown for a 
particular range of DM mass, the analysis can, in principle, be valid for any choice of DM 
mass. Since effective coupling determines the fundamental energy scale, A of the interaction 
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{Gi ~ ^), we can conclude from our study that the constrained value of such energy scale 
reduces as the mass of DM fermion becomes higher. 

In this work we thus have derived the self interaction parameters and their interrelations 
in an effective model for 4-point interaction among fermion WIMPs and further constrained 
the effective DM self couplings for cluster observational data. The derived constraints on self 
couplings Gi from self interaction are found out to be less severe by few orders of magnitude 
than those for the effective couplings among fermionic DM and SM fermions consistent with 
DM relic abundance bounds by Planck observation. The general analysis performed in this 
study can also be extended for other DM self interaction limits from other astrophysical 
sites such as from dwarf galaxies, Milky Way where the characteristic DM velocity scale is 
smaller than that of clusters. Alternatively since the DM particles are captured in the cores 
of massive astrophysical bodies like solar core, neutron stars etc. with capture rate prevailed 
by DM self interaction strength 99|, this can also be used to constrain DM self couplings in 
this framework. 

Although we have only studied the cases where each type of such interactions is the 
only effective interaction playing among the DM particles, we could have chosen the various 
combinations of all such interactions. This would constrain the the effective self couplings 
differently which can, in principle, be determined by combining the limits obtained in this 
work. Also there are scenarios with more than one WIMP on which effect of such self 
interaction constraints can be studied. Moreover the analysis performed for Dirac fermions 
in this work can be extended for Majorana cases also. In addition, the yielded constraints 
in this study from self annihilation will restrict several particle physics models for fermionic 
DM such as supersymmetric models with neutralino DM. 
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